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" ^ Abstract. Various sequences that possess explicit analytic expressions can 

04 be analysed asymptotically through integral representations due to Lindelof, 

^ which belong to an attractive but somewhat neglected chapter of complex 

^ analysis. One of the outcomes of such analyses concerns the non-existence of 

linear recurrences with polynomial coefficients annihilating these sequences, 
and, accordingly, the non-existence of linear differential equations with poly- 
nomial coefficients annihilating their generating functions. In particular, the 
corresponding generating functions are transcendental. Asymptotic estimates 
of certain finite difference sequences come out as a byproduct of the Lindelof 
I I approach. 

o 
u 

d Introduction 

S There has been recently a surge of interest in methods for proving that certain 

sequences coming from analysis or combinatorics are non-holonomic. Recall that 
a sequence (/„) is holonomic, or P-recursive, if it satisfies a linear recurrence with 

K*" coefficients that are polynomial (equivalently, rational) in the index n; that is, 

in d 

0^ ^PfcH/„_fc = 0, pk{n)eC[n], po^O. 

^ I k=0 

Put otherwise, its generating function f{z), called holonomic or D-finite, satisfies 
a linear differential equation with coefficients that are polynomial (equivalently, 

rational) in the variable z; that is, with /(z) := X]n>o /"-^"i 

k> ^(7fc(z)-^/(z) = 0, qk{z)eC[z], qe^O. 

Within combinatorics, the holonomic framework has been largely developed by 
Stanley, Zeilberger, Lipshitz, and Gessel P!l] HSJ [221 133 , who provided a rich 
set of closure properties satisfied by the holonomic class. Since the class of holo- 
nomic functions contains all algebraic functions, establishing that a sequence is 
non-holonomic can in particular be regarded as a strong transcendence result for 
its generating function. 
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In recent years, proofs have appeared of the non-holonomic character of sequences 
such as 

logn, ^fn, n", m^, log log n, C,(n), 

tin 



!, arctan(n), \/n^~+T, 



e 



where is a harmonic number, Wn represents the nth prime number, and C(s) 
is the Riemann zeta function. The known proofs range from elementary |16j to 
algebraic and analytic [21 [TOl |22l [23] . The present paper belongs to the category of 
complex-analytic approaches and it is, to a large extent, a sequel to the paper |10j . 

Keeping in mind that a univariate holonomic function can have only finitely 
many singularities, we enunciate the following general principle. 

Holonomicity criterion. The shape of the asymptotic expansion of a 
holonomic function at a singularity zq is strongly constrained, as it can 
only involve, in sectors of C, (finite) linear combinations of "elements" 
of the form 



(1) exp(p(Z-l/'■))z"^Q,(IogZ)Z^■^ Z:={z- 



zo 



j=0 



with P a polynomial, r an integer, a a complex number, s a rational 
of Q>o and the Qj a family of polynomials of uniformly bounded degree. 
(For an expansion at infinity, change Z to Z :— 1/z.) Therefore, any 
function whose asymptotic structure at a singularity (possibly infinity) 
is incompatible with elements of the form ([T]| must be non-holonomic. 

Equation ([ij is a paraphrase of the classical structure theorem for solutions of linear 
differential equations with meromorphic coefficients [33]; see also Theorem 2 
of |10| and the surrounding comments. 

What the three of us did in [10] amounts to implementing the principle above, 
in combination with a basic Abelian theorem. Functional expansions departing 
from can then be generated, by means of such a theorem, from corresponding 
terms in the asymptotic expansions of sequences: for instance, quantities such as 

loglogn, yi^, e^^, ... 

logn 

constitute forbidden "elements" in expansions of holonomic sequences — hence their 
presence immediately betrays a non-holonomic sequence. In proving the non- 
holonomic character of the sequences (logn) and {\/n), we could then simply 
observe in [TU] that the nth order differences (this is a holonomicity-preserving 
transformation) involve log log n and 1 / log n in an essential way. 

What we do now, is to push the method further, but in another direction, namely, 
that of Lindeldf representations of generating functions. Namely, for a suitable 
"coefficient function" (f>{s), one has 

1 <.l/2-|-ioo 

(2) J2^^n){-zr = -— <j>{s)z^-. ds, 

' ^ /ITT . Sin TT.Q 



2i7r 

n=l 



1/2-ioo 



where the left side is, up to an alternating sign, the generating function of the 
sequence {(j){n)). Based on representations of type we determine directly the 
asymptotic behaviour at infinity of the generating functions of several sequences 
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given in closed form, and detect cases that contradict (jlj, hence entail the non- 
holonomicity of the sequence (0(n)). 

Here are typical results that can be obtained by the methods we develop. 

Theorem 1. (i) The sequences e™ , with c, G M, are non-holonomic, except in 
the trivial cases c—Q or 9 £ {0, 1}. In particular, 

(3) e^, e-v^, e^/", e"!/" 

are non-holonomic. 
(ii) The sequences 

(4) t;^. -tVT^ nnV2), -P^, r(ra), 



2«±1' n! + l' ' ^(n^/3)' ' " C{n + 2) 

are also non-holonomic. 

The non-holonomicity of the sequences in ^ also appears in the article by Bell et 
al. [2], where it is deduced from an elegant argument involving Carlson's Theorem 
combined with the observation that (j){s) = e'^* is non-analytic and non-polar 
at 0. The results of P] do usually not, however, give access to the cases where the 
function (/)(s) is either meromorphic throughout C or entire. In particular, they do 
not seem to yield the non-holonomic character of the sequences listed in Q , except 
for the first one. (Indeed, 1/(2" ± 1) has been dealt with in pi by an extension 
of the basic method of that paper. Moreover, any potential holonomic recurrence 
for 1/(2" ± 1) can be refuted by an elementary limit argument, communicated by 
Frederic Chyzak and transcribed in [T71 Proposition 1.2.2].) 

A great advantage of the Lindelof approach is that it leads to precise asymptotic 
expansions that are of independent interest. For instance, in §2.2| below we will 
encounter the expansion 



^n! + l f-" smnskT'{sk + l)' 

n>l k>l 



OO, 



where 

isk)k>i ~ (-3.457, -3.747, -5.039, -5.991, . . . ) 
are the solutions of r(s -I- 1) = — 1. This formula features a remarkable structural 
difference to the superficially resembling function X]n>i(~'^)"/"' — exp(— z) — 1. 
Plan of the paper. The general context of Lindelof representations is introduced 
in Section [l] in particular. Theorem [2] of that section provides detailed conditions 
granting us the validity of ([2|. As we show next, these representations make it 
possible to analyse the behaviour of generating functions towards -|-c», knowing 
growth and singularity properties of the coefficient function (f)(s) in the complex 
plane. The global picture is a general correspondence of the form: 

location Sq of singularity of (j){s) — > singular exponent of F{z) 

nature of singularity of 4i{s) — > logarithmic singular elements in F{z) . 
(This is, in a way, a dual situation to singularity analysis [TTJ[T3].) Precisely, three 
major cases are studied here; see Figure [2] below for a telegraphic summary. 

— Polar singularities. When (j){s) can be extended into a meromorphic func- 
tion, the asymptotic expansion of the generating function of the sequence 
(^(n)) at infinity can be read off from the poles of 4'{s). Section |2] gives 
a detailed account of the corresponding "dictionary" , which is in line with 
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Coefficients (j){n) 



1 



2y¥(logz)i/4 l + z 

-^(logz)-i/''cos(2VI^-|^) ^ 



(1 + ^)3/2 ''''PV4(l + z) 



-1 + ^^ E(l)+E'(l)(l + z) 

VTT log Z 



Figure 1. Asymptotic forms of E(2;c, 6), for representative parameter values. 



early studies by Ford [TS]. It implies the non-holonomicity of sequences 
such as 1/(2" - 1), l/(n! + 1), T{n^/2). 

— Algebraic singularities. In this case, a singularity of exponent —A in the 
function (f){s) essentially induces a term of the form (log in the gener- 
ating function as z ^ +oo. We show the phenomena at stake by perform- 
ing a detailed asymptotic study of the generating functions of sequences 
such as e^ in Section [sj based on the use of Hankel contours. The non- 
holonomic character of the sequences e*" for 6 E ]0, 1[ arises as a conse- 
quence. 

— Essential singularities. The case of an essential singularity is illustrated 
by 4'{s) — e*^/*: in Section [Zj we work out the asymptotic form of the 
generating function at infinity, based on the saddle-point method. In this 
way, we also obtain the non-holonomic character of sequences such as e^^^^ 
by methods that constitute an alternative to those of 

As the discussion above suggests, the present article can also serve as a synthetic 
presentation of the use of Lindeldf integrals in the asymptotic analysis of generating 
functions. The scope is wide as it concerns a large number of generating functions 
whose coefficients obey an "analytic law". This is a subject, which, to the best 
of our knowledge, has not been treated systematically in recent decades (Ford's 
monograph was published in 1936). In particular, the joint use of Lindelof repre- 
sentations and of saddle points in Section [4] as well as the corresponding estimates 
relative to the family of functions 

oo 

(5) E(z;c,0) :=^e^"''(-z)", 

Tl=l 

appear to be new. Figure [l] summarizes some special cases of the results we obtain 
for E(z;c, 6*). Note that six essentially different expansions occur, depending on 
the parameter values and on the singularity we are interested in, and observe the 
surprising occurrence of subtle oscillations associated with e~^/". 

In Section [5] we show that the technology we have developed also provides non- 
trivial estimates in the calculus of finite differences. Finally, Section [6] completes 
our investigation of the function E(z; c, 9), by working out its asymptotic behaviour 
near z — —I. 
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1. Lindelof Representations 

Lindelof integrals provide a means to express a function, knowing an explicit 
law" for its Taylor coefficients. Let s i~+ 0(s) be a complex function that is an- 
alytic at all points of M>o; the (ordinary) generating function of the sequence of 
values {(l){n)) will be taken here in its alternating form: 

(6) F(z):=5]0(n)(-zr. 

n>l 

The function (j>{s), which is typically given by an explicit expression, represents the 
"law" of the coefficients of F(z): it extrapolates the integer-indexed sequence (0(n)) 
to a domain of the complex numbers that must contain the half-line K>i. The key 
idea is to introduce the Lindelof integral 

(7) K{z-C) -.^ ^ U{s)z 



2\iT Jq sin(7rs) ' 

where C is a contour enclosing the points 1,2,3, .. . and lying within the domain 
of analyticity of 0(s). Formally, as well as analytically (see Theorem [2] below) , 
when (/)(s) is well-behaved near the positive real line, a basic residue evaluation 
shows that, with e(C) = ±1 representing the orientation of C, 

(8) A(z;C) = e(C)^0(n)(-z)", 

n>l 

since the residue of tt/ sinvrs at s = n equals (—!)"■ Thus, the Lindelof integral ([t]) 
provides a representation of the generating function F{z) of (pi. Then, suitable 
growth conditions on (j){s) enable us to preserve the validity of (8| under suitable 
deformations of the contour C. We state: 

Theorem 2 (Lindelof integral representation). Let (j){s) he a function analytic 
in 5R(s) > 0, satisfying the growth condition 

(Growth) \(t){s)\ < C ■ e'^l^'l as \s\ ^ oo for some A e ]0,7r[ and C > 0, 

in 3fi(s) > 1/2. Then the generating function F(z) — X]ii>i '^("')(~-^)" ^•^ analyti- 
cally continuable to the sector — (tt — A) < arg(z) < (n — A), where it admits the 
Lindelof representation: 

Proof. By the growth condition, we have |0(«)| — 0(e"^"), so that F{z) is a priori 
analytic in the open disc \z\ < e~^. The proof proceeds in three moves. 

(i) Fix z to be positive real and satisfying z < e~^. Define the (positively 
oriented) rectangle TZ[m,N], with m,N e Z>o by its opposite corners at 1/2 — Ni 
and m -I- 1/2 -I- TVi. With the notation ([t]), Cauchy's residue theorem provides 

m 

A(z;7^K7V]) = ^0(n)(-z)^ 

n=l 

We first dispose of the two horizontal sides of this rectangle. For s in the complex 
plane punctured by small discs of fixed radius centred at the integers, one has 

(10) =ofe--l^^(^)l 
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A consequence of this estimate is that the integrand in the Lindelof representation 
decays exponentially with A^: for fixed z e ]0,e~'^[, we have 

0(s)z^^^ = O(e^^e"-^), 
sm TTS 

so that the integral along the two horizontal sides of 7?.[m, N\ is vanishingly small. 
One can accordingly let N tend to +00, which gives in the limit 

m 

A(z;7^[m,oo]) = ^0(r^)(-z)^ 

meaning that a partial sum of Fi^z) is expressed as the difference of two integrals 
along vertical lines. 

(ii) We next let m tend to infinity. With z still a fixed positive quantity satisfying 
z — e^^ for some B > A and s = m + 1/2 + it, we have, for some constants K, K' 

TT 



0(s)z" < ATexp AJ{m+ll2Y + 

sm vrs L 

< i4:'e('^~-^)"e(^"'')l*l. 



^-Bm^-K\t\ 



Thus, as a function of m, the Lindelof integral is 0(e'^'^^^''™). This shows that the 
contribution to the integral (|7| arising from the rightmost vertical side of 7?.[m, +00] 
is vanishingly small. Letting m tend to infinity then yields the representation Q 
in the limit (upon taking into account the change of sign due to orientation). 

{iii) Finally, the convergence of the integral persists, for any real z > 0, given 
the growth condition on 0: this ensures that F{z) is analytic at all points of the 
positive real line. Furthermore, for z = re"' and s — 1/2 + it, we have 

(^11-^ \-,s\ _ ^/^+'* 



so that the Lindelof integral ([9| remains convergent in the stated sector, where it 
provides the analytic continuation of F{z). □ 

This theorem was familiar to analysts about a century ago: it forms the ba- 
sis of Chapter V of Lindelof 's treatise |24j dedicated to "prolongement analytique 
des series de Taylor" and published in 1905; it underlies several chapters of W. 
B. Ford's monograph [15] relative to "T/ie asymptotic developments of functions 
defined by Maclaurin series''^ first published in 1936. Lindelof representations are 
also central in several works of Wright [HI about generalizations of the expo- 
nential and Bessel functions. Last but not least, this circle of ideas can also provide 
a basis to Ramanujan's "Master Theorem" [31 pp. 298-323], as brilliantly revealed 
by Hardy in [19l Ch. XI]. (We propose to return to properties of the associated 
magic duality''^ in another study.) 

2. Sequences with Polar Singularities 

The original purpose of Lindelof representations was to provide for analytic con- 
tinuation properties. For instance, as a consequence of Theorem |2] generalized 
polylogarithms, such as 



Lii/2,o(2) = ^ Lio4(z) = ^log) 
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are continuablc into functions analytic in the complex plane slit along the ray from 1 
to +00; see for instance [nHH]. Another fruitful corollary of the representations is 
the possibility of obtaining asymptotic expansions. In this section, we examine the 
simple case of generating functions whose coefficients admit a meromorphic lifting 
to C. 

2.1. Polar singularities. The following lemmgj^is used throughout Ford's mono- 
graph [15, Ch. 1]. 

Lemma 1 (Ford's Lemma, polar case). Assume that 0(s) satisfies the conditions 
of Theorem^ Assume that it is meromorphic in 5R(s) > —B and analytic at all 
points o/5R(s) = —B. Assume finally that the growth condition (Growth) extends 
to the larger half-plane 3fJ(s) > —B. Then, the generating function F(z) admits, as 
z — > +00, an asymptotic expansion of the form 

(12) F(z)^~ V Res( ^^(l3(s)z'';s = so) +0(z-^) , z ^ +00, 
^-^ \smns I ^ ' 

l/2>5R(so)>--B 

where Res is the residue operator and the sum comprises all poles of (f){s) / simrs 
that lie in the strip —B < 3fJ(s) < 1/2. 

Proof. Start from (lol, push the line of integration to the left, and take residues into 



account. This gives directly the expansion (12 1. □ 



The following observations are to be made concerning (12 1. 

(i) A pole of (/)(s)/ sin 7r(s) at sq and of order /i > 1 gives rise to a residue 
which is the product of a monomial in z and a polynomial in logz: 

z^" P{log z), where deg(P) = /i — 1. 

Such poles may arise either from 0(s) or from 1 / sin tts at s = 0, —1, —2, • • • . 
In the case where cjj^s) has no pole at s = —n, the induced residue is of the 
form (j){~n)z^^^ . 

(a) Additional poles of <j){s) in the right half-plane can be covered by an easy 
extension of the lemma, as long as they have bounded real parts and are 
not located at integers. 
(Hi) As a consequence of Item (i), poles farthest on the right contribute the 
dominant terms in the asymptotic expansion of F(z) at +00. 



(iv) The asymptotic expansions of type ( 12 1 hold in a sector containing the 



positive real line. (To see this, note that the Lindelof representation re- 
mains valid for z in such a sector and that the growth condition in an ex- 
tended half-plane guarantees the validity of the residue computation leading 



to (121.) 



2.2. Non-Holonomicity Resulting from Polar Singularities. We list now 
a few sequences that may be proved non-holonomic by means of Ford's lemma 
(Lemma [ij in conjunction with the non-holonomicity criterion based on ([l| . We 
restrict ourselves to prototypes; a large number of variations are clearly possible. 



^This lemma is of course closely related to its specialization to hypergeometric functions, of 
which great use had been made in early works of Barnes and Mellin; see |34| and |40l Ch. XIV]. 
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Point So <j}{s) F{z) 

Regular point 

So = — n, (l){s) analytic 4>(—''^) H (— 



So ^ Z 



Polar singularity [fj2j 

1 TT 

So ^ Z, simple pole of (h(s) z'^° 

S — So Sm TTSo 

So pole of order /i of 4'{s)/ sinTrs 2;''''P^_i(log2:) 
Algebraic singularity [^sj 

1 _ TT 2 = (log Z)^-^ 

(s — so)*^ sin ttsq r(A) 



Essential singularity [f|4j 

SoGC, 6I<0 ^Hs-so)" _A'i2=0(log^-)2(I^ exp(_R'2(log2)5^) 

soGC, 6»<0 e-(""''o^'' Liz"«(logz)5a^exp((L2 + iL3)(log2)5^) 



Figure 2 . Sample cases of the correspondence between local (regular 
or singular) elements of a function at a point so and the main asymptotic 
term in the expansion of the generating function F{z) at infinity. 

Proposition 1. The following sequences are non-holonomic (with i = ): 

I /l.n- T^, /2,„-r(nV2), h, 

(13) l + 

Proof, (i) First the sequences /i,n, /2.n, /s.n are treated as direct consequences of 
Lemma [T] 

For /i.„ = 1/(1 + nl), we observe that the extrapolating function (f){s) ~ 1/(1 + 
r(s)) is meromorphic in the whole of C. Thus the basic argument of 2, Th. 7] 
is not applicable. However, examination of the roots of the equation T{s) = —1 
reveals that there are roots near 



_ r(n%/2) 
^ C(n + 2)- 



-2.457024, -2.747682, -4.039361, -4.991544, -6.001385, -6.999801, -8.000024, -8.999997, 

and so on, in a way that precludes the possibility of these roots to be accommodated 
into a finite number of arithmetic progressions. (To see this, note that if T(s) = — 1, 
then r(l — s) = — Tr/sinTTs by Euler's reflection formula. As s moves farther to the 
left, the quantity r(l — s) becomes very large; thus sinTrs must be extremely close 
to zero. Hence s itself must differ from an integer — by a very small quantity, 
which is found to be ~ (— 1)'^^^/A;!, in accordance with the numerical data above.) 
When transposed to the asymptotic expansion of F{z) at infinity by means of 
Ford's Lemma (Lemma [T]), this can be recognized to contradict the holonomicity 
criterion ([T|). Indeed, the exponents of Z = z — zq that can occur in the singular 
expansion of a holonomic function are invariably to be found amongst a finite union 
of arithmetic progressions, each of whose common difference must be a rational 
number. 
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A similar reasoning applies to /2,„ = r(n\/2), /s.n = r(?i\/2)/r(n\/3), and 
more generalljj^ r(an), where a G M \ Q. For instance, in the case of /2,„, we 
should consider 0(s) — T{s^/2)/T{s)^ , where the normalization by T{s)^ ensures 
both the analyticity at of F(z) and the required growth conditions at infinity 
of Lemma [l] The poles of <j){s) are now nicely aligned horizontally in a single 
arithmetic progression, but their common difference (1/V2) is an irrational number. 

(m) Next, for the sequences /4,n, /s.n, /e.n, we can recycle the proof technique 
of Lemma [l] so as to allow for an infinity of poles in a fixed- width vertical strip 
(details omitted). 

In the case of /4^„, we are dealing with the function (j){s) = 1/(2^ — 1), which is 
meromorphic in C, and has infinitely many regularly spaced poles at s = 2ifc7r/ log 2, 
with fc € Z. The "dictionary" suggested by Figure [2] and Lemma [T| applies to the 
effect that F{z) has an expansion involving infinitely many elements of the form 
^2ife7r/iog2^ A similar argument applies to /s^n ~ r(in) which now has a half line of 
regularly spaced, vertically aligned, poles on 3?(s) = 0. 

As another consequence (but not a surprise!), the sequence /6,„ — 1/C("- + 2) is 
non-holonomic, since (j){s) — l/C(s + 2) satisfies the growth conditions of Theorem[2] 
and the Riemann zeta function has infinitely many non-trivial zeros. Likewise, 
for /4,n and /5.„, the expansion of F(z) contains exponents with infinitely many 
distinct imaginary parts. (Non-holonomicity does not depend on the Riemann 
hypothesis.) □ 

In summary, assuming meromorphicity of the coefficient function (j){s) in C ac- 
companied by suitable growth conditions in half-planes, sequences of the form ((j)(n)) 
are bound to be non-holonomic as soon as the set of of poles of (j){s) is not included 
in a finite union of arithmetic progressions with rational common differences. In a 
way this extends the results of [U |S] to cases of functions that are not entire. 



3. Sequences with Algebraic Singularities 

3.1. Analogue of Ford's Lemma. Ford also discusses the case where (j){s) has al- 
gebraic singularities. In this situation, it is no longer possible to move the contour of 
integration past singularities; Hankel contour integrals replace the residues of ( |12[ ), 
and only expansions in descending powers of log z (rather than z) are obtained |15l 
Ch. III]. See Figure |2] for an apergu. We first define the kind of singularities that 
can be handled. 

Definition 1. The function (j){s) is said to have a singularity of algebraic type 
(A, 6, tp) at So if a local expansion of the form 

(14) 0(s) = (s-so)-V((s-so)^) 
holds in a slit neighborhood of Sq, where A € C, ^{9) > 0, and 

(15) V(s) = EPfc«' 

fe>0 

is analytic at zero. 



The corresponding generating functions are related to classical Mittag-Lefiler and Wright 
functions |41l 1421 143| . They arise for instance in fractional evolution equations 28 and in the 
stable laws of probability theory [8| §XVII.6]. 
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To state the full result, we denote by 6j(so) the coefficients in the expansion 

(16) - V b,{so){s - so)^ so e C. 
sm TTs ^-^ 

j>-i 

Note that 6_i(so) = for sq ^ Z, while for rt G Z, we have 

(17) 6_i(n) = (-l)", 62fc(n) = 0, 62fc-i(n) = (-l)"(2-22-2'=)C(2fc). 

Lemma 2 (Ford's Lemma, algebraic case). Suppose that 4>{s) is analytic through- 
out C, except for finitely many singularities of algebraic type {Xi,9i,i})i) at Si, 
i — 1,...,M. The m-th branch cut should be at an angle uJm £ ] — ^tt, 0[ U ]0, |7r[ 
with the negative real axis, and the cuts may not intersect each other or the set Z>o. 
We impose the growth condition (Growth) with 

A < n min I sin LUrn I • 

l<m<M 

Assume furthermore that the singularities Sm dfe sorted so that 3?(si) = • • • = 
3?(sAr) > 5R(sAr+i) > • • • , for some I < N < M. Then F{z) has the asymptotic 



expansion 

(18) F{z)^ (-l)"+V(-n)z-" 

0<ri<Tlniax 

E..s^ Pm,kbj{Sm) ■,^-9,„fc-j + A„-l 

2^ r(-0„/fc-j + A„,)^^ ^ 

m=l k>0 ^ ™ J > mj 

where the Pm,k are the coefficients in the expansion of at Sm, in the sense 
of (15), and the summation range of the first sum in (18 I is determined by 



(19) 



-5R(si) 3?(si) e Z and s,„ (^Z for I <m < N 

[— 5R(si)] — 1 otherwise. 



The variable z may tend to infinity in any sector with vertex at zero that avoids the 
negative real axis. 

Proof. The statement assembles and generalizes results by Barnes [1] and Ford [15], 
who treat the case 6',„ = 1 (the latter reference offers a very detailed discussion). 
Extending the integration contour as usual (see Figure [s]) , we find the expansion 

(20) F{z) + — Y. H^)z'^^ds ^ (-ir+V(-n)z-", 

-n(^{si,...,SM} 

where 7i„ is a narrow Hankel-type contour, positively oriented and embracing the 
m-th branch cut. Suppose ffi'st that M = 1, si = 0, set bj :— bj{0), pk '.— pi,k 
(see (16l and (15l), and drop the index of etc. In a slit neighborhood of 

zero, we then have 

<^(.s)— = s-^ 5]p,6,s^'=+^-. 

sm TTS ^-^ 

k>0 

We first deal with a truncation 

TKis) s-^ Y P'^b,s'^+^ 

k>0,j>-l 
9k+j<K 
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Figure 3. A rectangular integration contour embracing branch 
cut singularities. 



of this expansion. Plugging Tr- in (20 1 and substituting y = — slogz, we then 
evaluate the integral termwise by Hankel's formula for the Gamma function: 



1 

'2i7r 



We thus obtain 
1 
2i7r 



Pkb] 



E 



\-yr+'-^dy 



fc>Oj>-l 
9k+j<K 



T{-ek - J + A) 



(logz)-H 

(logz)-«'=-^+^-i. 



We now have to show that 



(21) 



H 



\ sm TTS / 



For s inside the circle of convergence of if), we can expand 0(s)7r/ sin tts, substitute 
y — — slogz, and estimate the series tail X]efc+j>is' by the triangle inequality. 

To treat the remaining part of the contour 7i, we pick a point w on either of 
the rectilinear portions of 7i, and lying inside the circle of convergence of -0. By 
taking Ti narrow enough, both rays of 7i will admit such a point with ^{w) < 0. 
Then the portion J^^ °° of the integral in (21 1 is smaller than 



gSR(tolog2) 



(22) 

^ ek+j<K 

where the boundedness of \<j){s)'K / siuTrsj < C results from our growth assumption 
on(j)[s) and the fact that sin TTS = 0(exp(— 7r|3(s)|)). The integral in (22| converges, 
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and 

g5R(Milog2:) _ 0(|z|^'^"'^), 

which is neghgiblc in the logarithmic scale of the problem, as di{w) < 0. 

This completes the proof in the special case M — 1, si — 0. The full result 



follows from the special case by performing the substitutions s i— > s + Sm in ( 20 1 



Deleting redundant terms in the resulting expansion then yields (18 1. Note that 



in the first case of (19 1 we must include the contribution of the pole at s = 5R(si) 



in (18), which gives rise to the summand n — — 3fi(si). Otherwise, we need to 



consider only the poles whose real part is larger than 5R(si). □ 

We make the following remarks concerning Lemma [2] 
(j) Clearly, the statement extends to functions 0(s) having both poles and 



algebraic singularities. In fact the expansion (18 1 essentially remains valid 
then, as the reciprocal of the Gamma function vanishes at the non-positive 
integers. For instance, if s,„ ^ Z is a simple pole, then only the summand 
k = j = of the inner sum remains, which is in line with Lemma [T] The 



dominating singularities si, . . . , sjv that enter the expansion ( 18 1 must then 
not only comprise the rightmost algebraic singularities, but also the poles 
whose real parts are equal to theirs or greater. 

(a) We disallow horizontal branch cuts in the lemma, in order to take advantage 
of the exponential decrease of tt/ sinTrs along vertical lines. If a horizontal 
cut is present, and (j){s) stays bounded near the cut, the result persists. 

(Hi) There is a slight error in the statement of [TSJ §5]: Ford assumes that his 
function P, which corresponds to our (/)(— s)7rs/ sin its, is bounded in a right 
half-plane. This is usually too restrictive, due to the poles of l/sin7rs, and 
is in fact not satisfied by the application in [T5J §6]. 

(iv) By putting (f>{s) — s^'^ in Lemma [2[ we recover the classical expansion of 
the poly logarithm function at infinity [T5J . 

3.2. Asymptotic Analysis of the Generalized Exponential E(z; c, 9) when 
E ]0, 1[. Recall that in this case (j>{s) = exp(cs^), so that it is a straightforward 
application of Lemma [2] with Af = 1, A = 0, and ?A(s) = e'^". In fact this example 
was our initial motivation to extend Ford's result to the case where 0^ 7^ 1. As 
for the branch cut of the function , we may put it at any direction allowed by 
Lemma [2] The resulting expansion is 

(23) E(.;c,e)^- 5: fc.rSr-,) (^"g^)""^"' 

fc>o ■ ^ 
i>-i 



with bj{0) as given in (17). In particular, for the coefficient sequences e*^ we 
obtain 

E(z;c = ±l,0= i) = -lT^^ + of— 

V^rlogz \(logz)-'/^ 

3.3. Non-Holonomicity Resulting from Algebraic Singularities. The esti- 
mates of Equation ( 23 1 , when compared with the holonomicity criterion ([T]) , im- 



mediately yield the non-holonomic character of simple sequences involving the ex- 
ponential function, such as e*^ in Theorem [l] More generally, we can state the 
following result. 
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Proposition 2. Suppose that 4>{s) satisfies the assumptions of Lemma\^ and has 



a non-polar singularity at Si with an expansion of the type (14 1. Then the se- 
quence {(t>{n))„>i is not holonomic. 



Proof. This follows readily from the expansion ( 18 1 and the holonomicity crite- 
rion (jlj. Without loss of generality, we assume that si has maximal real part 
among the non-polar singularities of (t>{s). Now choose fco ^ such that pk^ 7^ 
and koOi — Ai ^ Z. (If there was no such fcoi then si would either be a pole or no 
singularity at all.) It then suffices to pick some jo > —1 with 6j„ ^ to exhibit a 
logarithmic term with "forbidden" , non-integral exponent and non-zero coefficient 



in(18l. □ 



Note that this proposition could also have been proved by the method of f2], 
based on Carlson's Theorem, but without the "constructive" feature of obtaining 
an asymptotic expansion of the associated generating functions. 

4. Sequences with Essential Singularities 

Going beyond the topics covered by Ford's treatise [15], we now investigate cases 
where the coefficient function (j){s) in ^ has essential singularities. We do not aim 
at a general statement here, but instead restrict attention to the generalized expo- 
nential E(z; c, 6) from ([s]), with 9 < 0. (In fact a general result encompassing both 
Theorems [3] and [4] below would probably be rather unwieldy.) This illustrates the 
use of Lindelof representations in conjunction with the saddle-point method [71ll3j. 
The resulting asymptotic formulas (Theorems [s] and |4] below) will be immediately 
recognized to be incompatible with the structure formula of ([T|): in this way, the 
present section completes our proof of Theorem[l] (Two rather easy supplementary 
arguments, which serve to cover the whole range of parameter values, but involve 



only crude asymptotic analysis, are collected in Subsection 4.3 ) 



4.1. Asymptotic Analysis of the Generalized Exponential E(z;c, 0) when 
c > and 6 < 0. In this section we determine the asymptotic behaviour of E(z) 
near infinity for positive c and negative 6. We present the analysis in the special 
case c = 1, 6 = —1. The generalization to arbitrary c > and 6* < is then easy. 
We start once more from the Lindelof integral representation 

1 <.l/2+ioo 

(24) Eiz;l,-l) = -— e^/'z'- ds. 



2i7r ,/1/2-ioo sinTTS 



Neglecting the effect of tt/ sIutts, the derivative 

as 

reveals a saddle point near s = where L := log \z\. We accordingly move the 



integration contour in (24) to the left, obtaining 



(25) E(z;l,-1) = -— / e^/'z'- ds. 

2i7r ' 



100 



We set s — L^^/^ + it. The main contribution to the integral arises near t — 0, 
say for \t\ < L^", where it will turn out that | < a < | is a good choice. The 
expansions which we require to approximate the integrand around the saddle point 
are collected in Figure |4] (We will recycle them in the next section.) From these 
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1 <"< ! 




exp(±i) = exp(±iLi/2 ± ^7,3/2^2 ^ .^^^^(i ^ 




- 1^1/2(1 + o(Li/2-")) 




sinTTS a V V 







Figure 4. Four elementary asymptotic expansions. The vari- 
able L tends to +00, and the first line specifies the range of s and 
the fixed parameters a, b, and a. 



we obtain, provided that a > |, the approximation 



= iLl/2g^p(2^1/2 _ ^3/2^2) . ^ 0(^2-3-)). 



This puts us in a position to evaluate the central part of the integral ( 25 1 : 

ds = / * dt 



2i JL-i/2_ii-Q sinTTS 27r 
Li/2e 



l/2„2Li/^ ./2L^/'*- 



1 i-3/4g-rV2d^ 



27r 7-y2L-V4-= ^ 

2L^''^ /-oo „ e^^^^^ 



(26) --^T^r^/ e-Z^dr^ 



2V27rLi/4 2y^Li/4 ' 

with a relative error of 0{L^^'^"). In order to let the integration bounds of the 
Gaussian integral tend to infinity, we have assumed a < | here. The tails of the 
Gaussian integral then decrease exponentially in L. 

In order to show that this is indeed the dominant part of the integral ( 25 1 , it 
remains to prove that the portion of the integral from ii^" to 100 (and thus, by 
symmetry, also from — ioo to — iL~") grows more slowly. 

First we consider t — 3(s) > 1. In this range we have e^/" — 0(1), and 

|z"| =exp(L^/2-targz) and l/siuTrs = 0(e-''*) 

lead to the bound exp(L^/^) • exp(— (tt + arg z)t)dt. To make the integral con- 
vergent, we assume that z tends to infinity in a sector that does not contain the 
negative real axis. 

Now consider L~" < t < 1. The factor is 0(exp(L^/^)) there, and the estimate 
|ei/-| < eVI-l ^ 0(exp(Li/2 _ ii3/2-2a)) 

follows from evaluating e^/'"' , which is a decreasing function of 3(s), at s = L^^/^ + 
i-L~". Since 

l/siuTTs = 0(l/s) = 0(L"), 
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we have established the tail estimate 



which grows slower than the absolute error in ( 26 1 . Hence the asymptotic behaviour 
of E(z; 1, —1) near infinity is 



z 00. 



2y7r(log2:) 

The error term follows from taking a G ] | , | [ close to | . 

All steps of the previous derivation are easily extended, when 1 /s is replaced 
by cs^, which yields the following result. 

Theorem 3. Let c > and 9 < be real numbers. Then 

(27) E(z; c, 9) = -Ki (log z) ^ exp(if2(log z) ) (l + O ((log z)"^) ) 

as z ^ 00 in an arbitrary sector with vertex at zero that does not contain the 
negative real axis. The positive constants Ki and K2 are defined by 

iCi (27r(l - 6l))~i/2(-ce')5(^fcr) and K2 '.^ {l - ^) {-c9)^ , 

and the exponent /i of the relative error estimate is 

-e 9>-2 



) 2(6-1) 



jhe 9<-2, 



with e an arbitrary positive real. 

4.2. Asymptotic Analysis of the Generalized Exponential E(z;c, 6*), when 
c < and 9 < 0. We present the detailed proof for the parameter values c ^ 9 = 
— 1. Then, the integrand of the Lindelof integral 

, .1/2+ioo 

(28) E(z;-1,-1) = -— / e-^/'^z'- ds 

2l7r7l/2-ioo SlUTTS 

has two saddle points, at ±iL^^/^, roughly, which will induce an oscillating factor. 
The argument of the axis of the upper saddle point is [7] 



7^ 1 d2 ^ 



A ' 



and that of the lower saddle point is jTt. We choose an integration path that has 
two segments passing through these saddle points at an angle of ±|7r with respect 
to the real axis and with length \/2L~°', where a is yet to be chosen. The segments 
are joined by a vertical line, and extended by vertical lines towards ±ioo. Our path 
thus consists of the five segments (cf . Figure [5| 

Ci : s = ~L-" +it, t< - L^", 

C2: s = -ii-'/2 + (l + i)t, N<i"", 

C3 : s = L-" + it, |t| < -i^", 

C4: s ^ iL^^^^ + {i - l)t, \t\<L-°', 

C5 : 3^-1-°+ it, t> L-^^^ + 1-°". 
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Figure 5. The landscape of |z''e^^/''/sin7rs|, where z = lO^*', and 
the new integration contour crossing the two approximate saddle 
points. 



We will see that the exponent a must satisfy the same bounds as in the previous 
subsection, i.e. | < a < |. For the segment C2, containing the lower saddle point, 
we again appeal to the expansions from Figure |4] and find 

- i7r-iLi/2 exp(-2iLi/2 - 2X^/2^2) . + QCi^-S")). 

sinTTS 

Here we have set s = —iL^^I'^ + (1 + i)i, and assume that a > |. Since 



exp(-2i3/2t2)d< ^ / e-'-'dr : 



for a < |, we thus have 



^^^ds = (i - l)(27r)-i/2L-i/V2'^''' • (1 + 0(L2-3")). 
sin ns 



The contribution of the upper saddle point, 

^'"""'^ -ds = (i + l)(27r)-i/2L-i/V'^''' • (1 + 0(L2-3o)), 



SlUTTS 



is similarly found. The dominant part of (28) is therefore 



2WC2UC4 sinTTs 2i V7c2 JC2J \JC2 sin'^s 

= -(27r)-i/2i-i/4(cos 2^1/2 + sin 2L'^^) • (1 + 0(L2-3")) 
(29) = -7r-i/2i-i/4 cos(2Li/2 - ivr) • (1 + 0(^2-3")). 

It remains to bound the integrals over Ci, C3, and C5. The portion of C5 with 
t = 3(s) > 1 is 0(exp(— L^"")), by the same argument as in the case of one saddle 
point. Now consider the lower part of C5, where we have L^^/"^ + L^" < < < 1. 
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The factor tt/ suitts is of order 0(i"). Also, it is easy to see that | exp(— l/s)| is a 
decreasing function of 3(s) there. At the lower endpoint of C5, we estimate 

exp(-l/s) = cxp{L^-" - 2L3/2-2a 

Since we have z" — 0(exp(— L^^")) in C5, this segment thus contributes only 
exp(— 2L^/2~2" +0(1)) to the integral. Finally, we examine the segment C3. The 
factor |exp(— l/s)| is an increasing function of |3(s)| there. Hence it suffices to 
estimate exp(— l/s)z* at the upper endpoint of C3, which is straightforward and 
shows that the integral over C3 is also negligible. This completes the tail estimate. 



Equation ( 29 1 hence yields the result 
E(z;-l,-l) = -^(logz)-i/4cos(2yb^-i^)+0((logz)-i/2+-), 
The generalization to arbitrary negative parameters is as follows. 
Theorem 4. Let c < and 9 < be real numbers. Then 

E(z; c,9) = Aiexp (^A2{log z)^^^ (logz)5(i^ cos (^As^logz)'^ + A^^ 

+ (exp (^A2(logz)^^ (logz) 50^"'"^ 

as z ^ 00 in an arbitrary sector with vertex at zero that does not contain the 
negative real axis. The constants are defined by 



A3 = (l-0-i)(c0)T^sin^, A4 = 

so that A2 is negative for —1 < 9 < Q, zero for 9 = —1, and positive for 9 < — 1. 
The exponent /i is as in Theorem^ 

Proof. The general proof is very similar to the special case c — 9 — —I (see above), 
upon taking into account the following comments. There might be more than two 
saddle points in general, but we have to consider only the ones that form a conju- 
gate pair having the largest real part, which are (approximately) at exp(±i7r/(l — 
5/))ii/(e-i)^ The saddle point axes have the arguments ±^7t{2 - 9)/{l - 9). When 
— 3 < < — 1, the proof proceeds as above. The parameter a, which governs the 
size of the two contour segments containing the saddle points, must satisfy 

<a< 2-« 



3(1-9) ^ " ^ 2(1-9) ■ 

For 9 < —3, there is one minor problem in the tail estimate: the decrease of 
(30) |exp(c/)| =exp(c|s|^cos(6larg(s))) 

does not hold on the whole of C5 . In fact it is plausible that there are oscillations 
here if \9\ is large. What happens is that, moving upwards along C5, the quantity 
I exp(cs^)| decreases up to a local minimum, at sq say. This results from elementary 
analysis; the argument of sq is arg(so) = 27r/(l — 9). Now the integral over the 
portion of C5 below sq can be estimated as before, by taking into account the length 
of this part of the contour and the size of the integrand at the lower endpoint of C5 . 
Above So, the crude estimate |exp(cs^)| < exp(— c|s|^) suffices. 

Finally, when —1 < 9 < 0, the contour has to be slightly adjusted, as the two 
saddle points at exp(±i7r/(l — 0))L^/(^~^) have negative real part, and we may not 
push the contour over the singularity at zero. Instead of joining the two central 
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segments C2 and C4 by a vertical line, we stretch each of them into the right half- 
plane, stopping at 5R(s) = L~°'. Then we join them by a vertical line, whose upper 
endpoint we call si. The value | exp(cs*)| decreases as we move on C4 to the right, 
until 3fi(s) = 0, from which point it increases. Moreover, as we move downwards 
from si, the quantity |exp(cs^)| decreases until 3(s) = 0. Therefore, it sufhces 
to check that the integrand at si is of growth slower than the central part of the 
integral, a property that is easily checked to hold. □ 

4.3. Completion of the Proof of Theorem [l} It remains to consider the pa- 
rameter region 9 > I, which is not covered by our previous asymptotic estimates 
of Fi{z;c,9). First, for parameter values c < 0, 1 < 6, which make E(2;;c, 6*) an 



entire function, a formula similar to (27 1 holds. The exponential growth order is 



the same as in (27), except for a different constant in place of K2. In particular, the 
sequence e™ is not holonomic for these parameter values either, as 6 /{9 — 1) 1. 
This asymptotic property is a special case of a result due to Valiron |38j . who used 
the Laplace method to investigate the behaviour of X]n>o as a; ^ 00, 

where G is a smooth function that satisfies certain regularity conditions. (Valiron's 
conditions actually require 1 < 6* < 2, but his analysis is easily extended.) 

Finally, in the remaining parameter range c > 0, > 1, the sequence e'^" grows 
faster than any power of n\, which is incompatible with the growth of any holonomic 
sequence. (In fact this observation shows that the (formal!) power series E(z; c, 9) 
does not even satisfy an algebraic differential equation |33j.) 



5. Asymptotics of Finite Differences 

Beyond establishing non-holonomicity, the analysis near infinity of generating 
functions, such as E(z; c, 9), is also of interest in the estimation of finite differences 
and related combinatorial sums. Given a sequence (/„), we shall refer to the derived 
sequence 

(31) Aj/]:-EQ(-l)Vfc 

as the sequences of differences. (In the standard terminology, we have Dn[f] = 
(— l)"A"/o; see [HlUni-) The relation between /„ and :— Dn[f] is translated 
at generating function level by the relation 

where f{z),g{z) are the "standard" (i.e., non-alternating) generating functions 
/(z):=^/„z", giz):=J29nz'\ 

n>0 n>0 

and F{z) = f{~z) — fo is the "alternating" generating function of (l6|. 



Because of the alternation of signs in (31 1 and the fact that the binomial coef 



ficients become almost as large as 2", the asymptotic estimation of differences is 
usually a non-trivial task. Here, the "surprise" is the fact that, for many explicit 
and simple sequences /„, the corresponding are much smaller than 2": huge 
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cancellations occur in (31 1 (see, e.g., [12], for cases related to /„ = n",logn, and 
so on). For instance, with /„ = e^^", n > 1, we find 

31 = -2.71828, 510 = -8.03246, <?ioo = -20.4159, giooo = -45.1379, 

and the sequence appears to grow rather slowly. For /„ — e^, it even appears 
numerically to tend slowly to 0. (For recent estimates relative to zeta values and 
inverse zeta values, see for instance, [14, .) We now explain how a Lindelof type of 
analysis can serve to quantify such phenomena. 

The basic message of singularity analysis theory jTTJ [13j |30] is that the behaviour 
of a sequence is (usually) detectable from the singularities of its generating function. 
Here, we should investigate the singularity of g(z) at 2 = 1. Now, this singularity is 
tightly coupled with the one of (1 — z)g(z), which by (32) depends on the behaviour 
of F(z) near +00. (This, by elementary properties of the conformal map z 1-^ w — 
z/(l — z) and its inverse w t-^ z = w/{l + w).) Clearly, by considering specific 
analytic maps cr(z) (here: (t(z) = z/(l — z)), a large number of seemingly hard 
alternating sums become asymptotically tractable. 

Corollary 1. The differences of the sequences e*^ and e^^^" have the following 
asymptotic behaviour. 



(33) 
(34) 
(35) 



E 

n 

E 

n 

E 



fc=i 



(_l)fce±^ 



(_l)'=e-i/fc 



±1 
-v/tt log n ' 

g2Vlogn 

'2V7F(logn)i/4' 
cos (2Vlog n - jTr) 



((logn)-i 



/4 



V^(logn)i/4 

Proof. Consider /„ = e*^, and define g{z) as above. Then (23) implies 

-1/2 



5(^) 



±1 1 



log 



1 



1 



1, 



whence (33 1 follows by the appropriate transfer theorem |TTJ [131 [30] . In the cases 



fn 



(with /o = 0), the growth of the difference generating functions. 



say gi{z) and g2{z), can be determined by Theorems |3] and |4] Note that both slowly 
varying and periodic functions can be subjected to singularity analysis fll[ll3[l37]. 
so that the formulas 




yield (34) and (35), respectively. 



□ 
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The estimate (34 1 bears a striking formal resemblance with the growth of the 
average value of the multiplicative partition function, which was found by Oppen- 
heim |31| and Szekeres and Turan [36j . Indeed, only the sign and the exponent 
of logn (— I instead of — |) differ. 

As an application of Corollary [T| we note that Madsen [27] has considered gen- 
eralized binomial distributions of the form 

n \ v-^ n ~ X 



where, for instance, he sets tt^ = exp((logp)fc°) with < a < 1 and < p < 1. We 



can then describe, by an obvious extension of (33 1, the way the probability mass 
function behaves for large parameters n: 



logp „ 
a; = 



r(l - a)(logn)'^ 
a logp 

a; > 1 



a;r(l - a)(logn)"+i 
In particular, there is no limit distribution, as n — > oo. 



6. Behavior of E(z; c, 6) at its Dominating Singularity 

Although not related to non-holonomicity or the Ford-Lindelof technique, it 
seems natural to complement the asymptotic results we have obtained for the func- 
tion ([5| by investigating its dominating singularity, located at z = — 1. The asymp- 
totic behaviour there is comparatively easy to determine. To begin with, for 9 < Q 
and any real c we can rewrite 

(36) E(.; c, 0) = ^ ^ ^ E ^ Li_,,(-z), \z\ < 1, 

n>lfc>0 ■ fc>0 

as a sum of polylogarithms 

n>l 

whose asymptotic behaviour at z = 1 is known. The shape of the asymptotic 
expansion of LIq, depends on whether a is an integer [9j [13] . 

Proposition 3. Let c he a real number and 6 be a negative real number. Then 
the asymptotic expansion of ¥i{z;c,6) at z = —1 is obtained by transporting the 
expansions oflA^ke into ( |36| ). 

Adding infinitely many asymptotic expansions termwise can be easily justified 
here by truncating the expansions and appealing to uniform convergence, which 
permits us to exchange limit and summation. For instance, if a > 1 is an integer, 
we have 

U^{z) = }^^w^-\\ogw~H^^,)+ V t3LQ^a-j)w\ 
(a — 1)! ^ V. 
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where w = — logz and Ha-i is a harmonic number. For the parameter values 
c=l,0=— Iwe thus obtain 



E(z; c = 1, = -1) ^ ^ ei/"(-^)" = E ^ Lifc(-z) 



n>l fc>0 

= + log ^ + C + 0((1 + z) log -^), 
1 + z \ + z I + z 

where C = -1 + Efc>2 C(fc)/fc! ~ 0.078189. 

We proceed to the case c > and < 6* < 1. For fixed z inside the unit disk, 
the summands have a peak at some n and then decrease rapidly, which makes 
the Laplace method a natural approach for estimating the sum. Our illustrative 
example is c = 1 and 9 = ^. The summands e^(— z)" have their peak near u :— 
l\v~'^\, where v :— — log|z|. This leads to the lower bound exp{jv~^){l + 0{v)), 
which was already noted by Borel ^ p. 69]. It is straightforward to determine the 
second order approximation of the summand near n — u, and to estimate the tails 
of the original sum and the second order approximation. What we find is 

E(z;c=l,0=i) = Eev^(-z)" 

n>l 

The (again straightforward) generalization from ^/n to cn^ reads as follows. 
Proposition 4. Let c be a positive real number, < 9 < 1, and e > 0. Then 

(37) E(z;c,e') = Ci(l + z)5I^exp(C2«i^)(l + 0((l + z)'^)) 

as z tends to —1^ in M, where v :— — log \z\. The constants Ci, C2, C3 are positive 
and given by 

Ci ^/2n{\-9)-^'^{c9)^^^^ , C2 := ^{c9)^ , 

C3 ic9)^, 

and the exponent of the relative error estimate is 

A* :=min{2(i^ -£,!}• 
Finally, we consider c < and < 6* < 1. Then the series 

E{z;c,9) = J2e^"\-zr 

n>l 

converges at z = —1. Inside the unit circle we may differentiate it termwise arbi- 
trarily many times, and the series thus obtained converge at z = —1, too. By Abel's 
convergence theorem, these values equal the limits of the derivatives as z — > — 1"*". 
The (divergent) formal Taylor series of E(z) at z — —1 obtained in this way is an 
asymptotic series for the function [TQ p. 30], which yields the following result. 

Proposition 5. Suppose that c < and < 9 < 1. Then 

(38) E(z;c,6') - uq + ui{l + z) + U2il + zf + . . . 
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as z tends to —1^ in M, where the coefficients are given by 
(39) Uk 



Note that the series (38 1 does not converge in any neighborhood oi z — — since 
Stirhng's formula yields 



where n(A:) = [(— fc/c0)^/^J approximates the index of the largest summand in (39 1. 
Hence z = — 1 is indeed a singularity. 

7. Conclusion 

We have revisited a classical method for the analytic continuation of power se- 
ries beyond their disc of convergence, with the goal of obtaining asymptotic expan- 
sions and comparing them against the possible expansions of holonomic functions. 
Our estimates can be used as building blocks for the asymptotic analysis of more 
complicated functions than those we have explicitly mentioned. For instance, the 
expansion of functions such as 

^ e^^(-z)" 

n>0 

at -|-c» readily springs from Lemma [2] in fjs] and series in the spirit of 

gl/n+l/ySj-__^-jn 

n>0 



can be analysed similarly as the ones in S4.1 and S4.2 As regards proving non- 
holonomicity, our results compete with those of Bell et al. 2^ , who also deal with 
sequences having an analytic lifting. Roughly speaking, our approach is more ver- 
satile for meromorphic functions, equivalent in the algebraic case, and less flexible 
in the presence of essential singularities. 

Note, however, that neither we nor Bell et al. [2] can show non-holonomicity 
of sequences whose extrapolating function is entire. For instance, we leave the 
non-holonomicity of sequences like cos(y^) and cosh(y^) as an open problem, 
since their analytic liftings, namely, cos(y^) and cosh(y's), have no singularity at 
a finite distance. 
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